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Abstract. We study the A'A' interactions in the liJ") = 0(1/2") and 1(1/2") channels and associated exotic 
state 0* from 2+1 flavor full lattice QCD simulation for relatively heavy quark mass corresponding to m„ = 
871 MeV. The s-wave KN potentials are obtained from the Bethe-Salpeter wave function by using the method 
recently developed by HAL QCD (Hadrons to Atomic nuclei from Lattice QCD) Collaboration. Potentials in 
both channels reveal short range repulsions: Strength of the repulsion is stronger in the 1=1 potential, which is 
consistent with the prediction of the Tomozawa- Weinberg term. The / = potential is found to have attractive 
well at mid range. From these potentials, the A^A' scattering phase shifts are calculated and compared with the 
experimental data. 



1 Introduction 

First evidence of the 6*^(1540) (an exotic state with baryon 
number B - I and strangeness 5 = -i-l which corresponds 
to the quark content uudds) has been reported by the LEPS 
Collaboration at SPring-8 1 1 1. Although numbers of exper- 
imental studies have been performed since then, the exis- 
tence of the 6>+(1540) is still controversial: CLAS Collab- 
oration observed no signal from their high statistic data f2l, 
and other experiments at high energies with high statistics 
and good mass resolutions did not find positive evidence 
as reviewed in 1 3 1 . On the other hand, LEPS Collaboration 
and DIANA Collaboration have recently reconfirmed the 
6*^(1540) signal with high statistics |4 5|. They indicate 
that the production mechanism of the (9^(1540) might be 
highly reaction dependent if it exits in nature. 

Theoretically, QCD studies of the mass and the quan- 
tum numbers of the 0^(1540) have been attempted using 
lattice QCD simulations and QCD sum rules. In the lat- 
tice QCD studies with the quenched approximation, ex- 
istence of the low-mass pentaquark is not conclusive yet 
(see e.g. the summary given in |6|). Recent QCD sum rule 
studies |7| with the operator product expansion up to di- 
mension 14 show some candidates of the 0^(1540) in the 
liJ") = 0(1/2-), 1(1/2-), 0(3/2+) and 1(3/2+) channels. 
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Determination of the quantum numbers of 6>+(1540) was 
also attempted from the phase shift analyses of the KN 
scattering data, and some candidates in L2/,2y = £^03,^05 
were reported ||8J. 



The main purpose of the present study is to investi- 
gate the low energy KN potentials in the UJ") - 0(1/2-) 
and 1(1/2 ) channels and to shed new light on 0+ from 
the 2+1 flavor full QCD simulations. A systematic method 
to extract the hadron-hadron potential from the equal-time 
Bethe-Salpter amplitude measured on the lattice has been 
recently developed and applied to the nucleon-nucleon po- 
tential by HAL QCD Collaboration 191101111121131 The 
potential obtained in this method can be utilized to calcu- 
late the scattering observables and to study the resonances 
and bound states. We utilize this method to extract the KN 
potentials and phase shifts. 



This paper is organized as follows. In section 2, the 
formalism to extract the KN potential from lattice QCD is 
briefly reviewed. Our numerical setup of the lattice QCD 
simulation is then shown in section 3, and the results are 
shown in section 4. These results are discussed in section 
5, and summary is given in section 6. 
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2 Formalism 



3 Numerical setup 



Following the basic formulation to extract the nucleon- 
nucleon interaction |9 10], we briefly show the equations 
to obtain the KN potentials below. We start with an effec- 
tive Schrodinger equation for the equal-time Bethe-Salpeter 
(BS) wave function 0(r): 



V2 r 



drU(r,r')cf>(r') ^ E4>(r), (1) 



where fi(- micmi^/(mic + niM)) and E denote the reduced 
mass of the KN system and the non-relativistic energy, 
respectively. The non-local potential U (r, r') can be ex- 
panded in powers of the relative velocity v = V/yU at low 
energies, 

t/(r, r') = y(r, \)6{r - r') 

= (Vwir) + VMLoir) + ■ ■ -Wr - r'), (2) 

where the N"LO term is of order 0(\"). In the leading or- 
der, we have 



V(r) - VLo(r) = —m + E. 



(3) 



In order to obtain the BS wave function of the KN sys- 
tem on the lattice, let us consider the four-point correlator: 

^,(x, y, t - to; J") = <0| Kix, t)NAy, OJKNito; J") |0) 
- 2 A„ (01 K(x, t)NJy, t) In) g-^-C-'o), (4) 



with the matrix elements 

A„^{n\JKN{to\J''m. 



(5) 



Here STKNito', J") denotes a source term which creates the 
KN system with spin-parity J" on the lattice. The four- 
point correlator in Eq. ^ is dominated by the lowest en- 
ergy state with total energy £0 at large time separation 
(f » fo): 

g„(r,t-kuJ'') = ^^„(x,y,f-fo;n 
y 

^ Ao0„(r; (6) 

with r = X - y. Thus, the KN BS wave function is de- 
fined by the spatial correlation of the four-point correlator. 
In Eq. (|6]l, we assume the Dirichlet boundary condition in 
temporal direction, so that the temporal correlation has an 
expnential form, g^^oC^'o). 

The BS wave function in s-wave state is obtained under 
the projection onto zero angular momentum (f '^'"°'), 

0(r; 1/2-) = 4 Z °V.(0"'r; 1/2-), (7) 

where g & O represent 24 elements of the cubic rotational 
group, and the summation is taken for all these elements. 
Using Eq. (O and Eq. dTji, we will find the KN potential 
and wave function from lattice QCD. 



In order to calculate the KN potentials in isospin 7 = and 
7=1 channels in 2h-1 flavor full QCD, we have utilized 
gauge configurations of JLDG(Japan Lattice Data Grid)/ 
ILDG(International Lattice Data Grid) generated by CP- 
PACS and JLQCD Collaborations on a 16^x32 lattice II14I15L 
The renormalization group improved Iwasaki gauge action 
and non-perturbatively 0(a) improved Wilson quark ac- 
tion are used at y6 = L83 which corresponds to the lat- 
tice spacing a = 0.1209 fm with the p meson mass in 
the chiral limit. The physical size of the lattice is about 
(2.0 fm)^ and the the hopping parameters are taken to be 
/f„ =Kd^ 0.1378 and k, = 0.1371. 

In the present simulation, we adopt the spatial wall 
source located at fo with the Dirichlet boundary condition 
at time slice f = fo -H 16 in the temporal direction and the 
periodic boundary condition in each spatial direction. The 
Coulomb gauge fixing is employed at f = fo. The number 
of gauge configurations used in the simulation is 700. 



4 Results 

The masses of hadrons are obtained by fitting correspond- 
ing two-point correlators. The obtained masses of the pion, 
kaon and nucleon are m;r = 870.7(1.9) Me V, = 911.5(1.9) 
MeV and mm - 1795.5(6.9) MeV, respectively. Thus the 
KN threshold energy is 2707 MeV. 
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Fig. 1. Effective mass plot of the KN states in / = O(red) and / = 
l(blue) channels. The dashed line denotes the threshold energy 
of the KN measured in this simulation. 



Fig. [T] shows the effective masses of the KN states in 
the 7 = and 7=1 channels together with the KN thresh- 
old energy. Eff'ective masses are obtained from the tem- 
poral correlation in Eq. (|6]l. We observe plateaus at large 
time separation, f - f > 7. The best fit in the plateau gives 
M/=o = 2708(1 1) MeV and M/=i = 2761(10) MeV, so that 
the lowest KN state is the 7 = channel. 

Figs.|2]and[3]show BS wave functions of the KN scat- 
terings in the 7 = and 7=1 channels, respectively. They 
are obtained from the lattice QCD simulations at large time 
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Fig. 2. The BS wave function of the KN scattering in the / = 
channel. 
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Fig. 4. The potential of the A^A' state in the / = channel without 
the energy shift E in Eq. 
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Fig. 3. The BS wave function of the KN scattering in the 7=1 
channel. 
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Fig. 5. The potential of the A^A' state in the 1=1 channel without 
the energy shift E in Eq. 



5 Discussion 



separation, f - f o = 8. The large r behavior of the BS wave 
functions in both channels do not show a sign of bound 
state, though more detailed analysis is needed with large 
volumes for a definite conclusion. The small r behavior of 
the BS wave functions suggests some repulsive interaction 
at short distance (r < 0.3 fm). Also, the repulsion in the 
7=1 channel seems to be stronger than that in the 7 = 
channel. 

The potential V{r) without the constant energy shift E 
in Eq. ([5]) for the 7 = (7 = 1) KN state is shown in 
Fig- El (Fig -El). These potentials are also calculated by us- 
ing the data at f - fo = 8. As expected from the BS wave 
functions in Figs.|2]and[3] we observe the repulsive interac- 
tions at short distance in both channels. Also, we observe 
the attractive well in the mid range (0.4 < r < 0.8 fm) 
in the 7 = channel. In the constituent quark model of 
hadrons [16|, similar short range repulsion in KN system 
has been predicted, while the attraction has not been found. 



Results shown in the previous section indicate that there 
are no bound states in the liJ") = 0(1/2") and 1(1/2") 
states for the pion mass m„ ~ 870 MeV. 
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Fig. 6. The potentials multiplied by a volume factor {p-V(r)). 
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In order to compare the strength of the repulsion be- 
tween the different isospin states, we plot the potentials 
multiplied by a volume factor (r^V(r)) in Fig. |6] with the 
energy shifts E = -5.0 MeV for the 7 = channel and 
£■ = 5.0 MeV for the 7=1 channel. These energy shifts are 
estimated from the asymptotic behavior of the potentials at 
large r, though analyses with large volumes are needed to 
determine these energy shifts more precisely. As seen in 
Fig. |6] the repulsion at short distance in the 7=1 chan- 
nel is much stronger than that in the 7 = channel. This 
is anticipated from the Tomozawa- Weinberg (TW) term in 
the effective chiral Lagrangian of mesons and nucleon: The 
7 = interaction vanishes, while the 7=1 interaction is 
repulsive from the contact TW interaction, which should 
be compared approximately with the integral of our r^V{r) 
at short distance. 
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Fig. 7. The phase shifts of the KN scatterings as the function of 
the momentum (Plab MeV) in the laboratory frame. The red solid 
(blue dashed) curve shows the phase shift of the I = 0(1 = 1) KN 
scattering. Data are taken from Ref. 1171 . 



By using the potentials which fit the lattice data in Figs.|4] 
and|5] we can calculate observables such as the scattering 
phase shifts. Fig.Qshows such phase shifts of the KN scat- 
tering together with the experimental data as a function of 
the laboratory momentum of the kaon. Theoretical phase 
shifts are calculated from the Schrodinger equation with 
the potentials, V{r), where the energy shift Ts = -5.0 (5.0) 
MeV in the 7 = (7 = 1) channel is taken into account. 
Although the hadron masses are heavy in the present sim- 
ulation, qualitative behabior of the phase shifts and also 
the relative magnitude between the 7 = and 7=1 chan- 
nels are consistent with the experimental data. Simulations 
along this line with ligher quark masses would eventu- 
ally lead to a definite conclusion on the existence of 0^ 
in li J") = 0(1/2-) and 1(1/2") channels. 



6 Summary 

We have performed the 2+1 flavor full lattice QCD simu- 
lation to investigate the KN interaction and a possible sig- 
nature of the exotic resonance 0^. The s-wave 7 = and 



7=1 potentials are extracted from the BS wave functions 
for relatively heavy quark mass corresponding to = 871 
MeV. Potentials in both channels reveal short range repul- 
sions: Strength of the repulsion is stronger in the 7=1 
potential, which is consistent with the prediction of the 
Tomozawa- Weinberg term. The 7 = potential is found 
to have attractive well at mid range. From these potentials, 
the KN scattering phase shifts are calculated and compared 
with the experimental data. Although the quark mass is 
heavy in the present simulation, the results indicate that 
our method is promising for future quantitative studies of 
the KN interactions and the exotic resonances in lighter 
quark mass region. 
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